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Abstract 

We calculate the electroweak and QCD corrections to W~W + — > tt and ZZ — > ti . We also 
consider the interplay of these corrections with the effect of anomalous interactions that affect the 
massive weak bosons and the top. The results at the VV level fusion are convoluted with the help of 
the effective vector boson approximation to give predictions for a high energy e + e~ collider. 
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1 Introduction 



The fact that the Higgs has not been discovered yet is evidence that the mechanism behind electroweak 
symmetry breaking, EWSB, giving masses not only to the weak vector bosons but also to fermions 
remains a mystery. Moreover the implementation of the Higgs within the Standard Model, SA4, poses a 
few conceptual problems. This has triggered an intense activity in model building in order to solve the 
hierarchy problem leading to a rich phenomenology at the upcoming high energy colliders, such as the 
Large Hadron Collider (LHC)[TJ and the International Linear Collider (ILC) [2]. 

In practically all the New Physics models which address the hierarchy problem the top plays a 
central and crucial role. One can mention models of dynamical electroweak symmetry breaking in the 
technicolour vein where the concept of top condensation [3] is exploited. In supersymmetry, the top keeps 
the the so-called minimal model alive through radiative corrections to the Higgs while in mSUGRA, 
for example, the top can trigger electroweak symmetry breaking. The phenomenology of the top in 
models of extra-dimensions is also special as discussed very recently [3]. It is no wonder that the top 
plays such a central role in probing the sector of EWSB as well as the flavour problem. After all, its 
large mass[5j is of the order of the spontaneous electroweak symmetry scale vj\f2 = 17AGeV which, 
in the SA4, gives it a Yukawa coupling of order one, incidentally sensibly of the same strength as the 
QCD coupling constant at the EWSB scale. This also explains why radiative corrections involving the 
top, especially in processes related to symmetry breaking, may be competitive with QCD corrections. 
It is in this context that we compute in this paper the electroweak corrections to vector boson fusion 
into a pair of top quarks, ZZ —* ti and W~W + — > ti. Weak vector boson fusion is the privileged 
process revealing spontaneous symmetry breaking through the longitudinal or Goldstone component of 
the vector bosons. Although most analyses have dealt with scattering into vector bosons, it is clear 
that fusion into tops is also a probe of symmetry breaking. Since the effects of the EWSB due to New 
Physics can be subtle, it is important that we know quite precisely the SA4 prediction. This calls for 
calculations beyond tree-level. QCD corrections have been performed recently in We will include 
here the electroweak corrections alongside the QCD corrections. We also compare the loop corrections 
effects with those from new operators parameterising the EWSB effects in W~W + — > ti. The complete 
study of vector boson fusion at the colliders would require a full 2 — > 4 calculation. Especially for 
the inclusion of radiative corrections this is a most daunting task. Though most of the machinery has 
been developed to tackle such complex 2 — > 4 processes at loop level, there have been till now only two 
such computations [8]. It therefore seems fit in a first step to revert to some approximation. The 
effective vector boson approximation, EVBA[9j, seems suited for such processes. This approximation has 
been extensively used in the context of weak vector fusion. We will conduct our calculation within this 
approximation for application to the e + e~ linear collider. The important genuine electroweak corrections 
should be captured in the subprocess. We only restrict ourselves to the linear collider; W~W + —* ti in 
the context of the LHC suffers from considerable background [10] mainly due to the overwhelming direct 
ti production. The Higgs range we consider in this paper does not go beyond a Higgs mass, Muiggs-, of 
300 GeV in accordance with the current limit from precision measurements [TT] . We intend to study the 
radiative corrections in the presence of a very heavy Higgs in a forthcoming publication. 

This paper is organized as follows. In the next section we briefly review ZZ — > ti and W~W + — > ti and 
show the importance of the longitudinal modes especially at high energies. We then expose the EVBA. In 
section 3 we present the main ingredients that enter the calculations at one-loop, the checks we performed 
as well as the tools we exploited. Section 4 gives our results for the electroweak and QCD corrections. 
Section 5 shows these results after convolution with the EVBA structure functions for the linear collider. 
Section 6 introduces the anomalous operators and their effects on the processes we study. Section 7 
summarises our conclusions. An appendix lists all formulae for the helicity amplitudes at tree-level, in a 
compact way, including the anomalous operators. 
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2 Tree-level overview and dominance of the longitudinal modes 
2.1 Notation and Born contribution 




(a) ZZ -> tt 





(b) W~W + -> tt 



FIG. 1: Tree level Feynman diagrams for the process a) ZZ — > tt and &j PF IF + — > tt. 
The tree level Feynman diagrams for the processes we study, Z{k\,\\) + Z{k2,^2) - » t(ks, 



A 3 ) + 

t(&4, A4) and W~(ki, Ai) + W + (k2, A2) — ► i(&3, ; A3) + t(fci, A4) are shown in figure 1(a) and 1(b) The 
labelling fcj and Aj = ±,0 stand for the momentum and the helicity. Complete helicity amplitudes 
formulae for W~W + — > tt are given in the Appendix. 

We are mostly interested in the longitudinal modes of the vector bosons, with Ai = A2 = = L. 
Total cross sections for this mode are shown in Figs EJ 
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FIG. 2: Total Born cross section for the processes W L — ► tt and Z^Zl — > tt for Muiggs 
1 20. 200, 300 GeV and p^ > 10 GeV. 

Z^Zl is more sensitive to the Higgs mass dependence especially for low centre of mass energies. The 
W mode provides though more yield. Note also, that for the range of Higgs masses we are studying the 
cross sections after reaching a maximum around threshold, drop as the energy increases. This should 
be taken into account when considering the effect of New Physics (NP) that tends to increase the cross 
section as the energy increases. When convoluting with the "structure function" of the EVBA (see next), 
the lower end of the cross section might pollute the NP effects. One should therefore aim at imposing 
a cut on the invariant mass of the system. We will always take a cut on the invariant mass of the tt, 
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m t f > 400GeV as is done in[6]. This cut will help improve the EVBA. We will also impose a cut on 
the transverse momentum of the top, > lOGeV. We will do this not only after convolution on the 
structure function, i.e. at the e + e~ level but also at the VV — > ti level. The figures we show are with 
these cuts. The p 1 ^ cut helps reduce some QCD background and fake pr in tij, see[T2J [THJ H3]. 
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FIG. 3: Relative rates for the Born cross section for different polarisations of the incoming vector bosons 
with respect to the case where the incoming bosons are both longitudinally polarised. For WW, AjAj 
labelling of the helicities refers to W~ (\{)W + . 



Another reason for aiming at a higher cut on the invariant VV (or similarly tt) system can be gleaned 
in Fig. [3J First VlVl scattering dominates for all energies. This said, at the lowest centre of mass energies, 
some transverse contributions might not be so negligible as compared to VlVl fusion, though they drop 
rather rapidly as the energy increases. Note however that the longitudinal/transverse (VlVt) mode in 
the SAi might bring some contamination at all energies, which could worsen the EVBA approximation. 



2.2 Born angular distributions 

Radiative corrections and the New physics can not only change the total yield but can also, and most 
dramatically, distort the shape of kinematical distributions. We show in Fig. 0] the angular distribution 
in W~W + — > tiior yjsww = 400 GeV, not too far from threshold and at high energy yjsww = ITeV. 
We choose to show W~W + — > ti because of the chiral structure of the W coupling and the small mass 
of the bottom (exchanged in the t-channel). This leads, at high energy, to an overwhelming production 
of top in the backward region, with respect to W~. In ZZ, the distribution are of course symmetric 
and much less peaked. A cut on the forward backward region will also help bring out the signal of New 
Physics which occurs at higher pr- 



3 



A 



M„ w = 1 20 GeV 
Vs_,=400 GeV 



Vs„= 1 000 GeV 



X. 



\ 



0.2 0.4 0.6 0.8 1 

cos(i>) 



-I -0.8 -0.6 



-0.2 0.2 0.4 0.6 0.8 1 

cos(tS) 



iggs 



120 GeV at 



FIG. 4: Angular distribution for the unpolarised W~W + — * tt cross section for Mh 
y/ s ww = 'iOOGeV and ITeV. 6 is the angle between the W~ and the t. The pjt > 10 GeV corresponds 
to a too small 9 cut that can not be represented on the plot. 



2.3 Weak bosons effective distribution functions 



To turn the subprocess ZZ — > tt and W~W + — ► tt into e + e~ — ► tie + e~ and e + e~ — > ttv e v e respectively, 
we revert to the EVBA approximation [9] . The approximation introduces the notion of the W/Z content 
of the electron (and positron). At high energies, we take the structure functions of the weak vector bosons 
inside a fermion as given by [151 116j : 



ffermion/V,\(.^ 



1 



16-7TX . 



(gv - ^9a? + (l - x) 2 ( gv + \g A f 



In 



p2 

T,max 



A = ±1 



(2.1) 



f fermion/V,\ (jE) 



A = 0, 



(2.2) 



4tt 2 x 

x is the momentum fraction transferred to the vector boson by the initial fermion. gy and gA are the 
vector and axial-vector couplings of the weak vector bosons to fermions. The transverse distribution is 
subject to more uncertainty due to the logarithmic factor. In our case, we take PT,max the maximum 
transverse momentum of the vector bosons allowed by the kinematics. In the collision of two elementary 

(3) (3) 

fermions f\ and fi of respective charge and isospin assignement \Qfc,T\ ) and {Qf 2 ,T^ ), the effective 
luminosity of the longitudinal vector bosons is given by: 

dC , , f 1 dx 



dr 



hh/v L v L 
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This gives for example for WlWl[17]: 

dC 

~d7 



(r) 



a 



47TS2, 



hh/w L w L 

The total cross section for the processes (e + e 



2 1 / 1 \ 
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^—^ ti + v e v e ) is then given by: 
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(2.5) 
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Formulae for the ZZ case can be derived in a similar manner. 

After convolution, the e + e~ cross section to ti through vector boson fusion is reduced much and 
will most probably be exploited only in a second stage linear collider with energy in excess of ITeV, see 
Fig. Note that at the level of e + e~ , the litte Higgs mass dependence there is in the WW channel 
(Mffiggg = 120GeV and Muiggs = 300GeV) is washed out but not in the ZZ channel, though the latter 
cross section is now, at least, an order of magnitude smaller than in WW suffering from the smaller 
couplings of the Z to the electrons. 




FIG. 5: Born total cross sections for (a) e + e — ► tte + e through Z^Zl — > tt and (b) e + e — > ttuu 
process through W£W^ — ► ti. 

One should keep in mind that the effective EVBA works best at high energy (compared to the W mass 
for instance). Reference|13j made a comparison of the EVBA with a full SM calculation (at tree-level) 
using CompHEP[18j. It was found that the EVBA is a good approximation for yjs ee = 1.5TeV and above. 

3 Electroweak and QCD corrections: 
Set up and details of the calculation 

3.1 General structure 

The calculation of the complete one- loop electroweak and QCD corrections to W~W + — ► ti and 
ZZ — > ti is performed with the help of two automatic codes for loop corrections, SloopS [19] and 
FormCalc/LoopTools[20j. Although SloopS uses many modules of [20J, the model files and the genera- 
tion of the Feynman rules is generated automatically with the help of LanHep[2"T]. Moreover, some loop 
integration routines have been improved. The strength of SloopS is the ability to perform very powerful 
gauge parameter independence checks as we will see later, a feature that is found in GRACE-loop[22j. 
SloopS has been developed for supersymmetry but has also a SM module. The results of SloopS, 
FormCalc/LoopTools and GRACE-loop have been checked for a variety of electroweak processes. All 
codes adopt the on-shell renormalisation scheme according to[22 } 1231 124j . In particular we take as input, 
the masses of the particles in the model and the electromagnetic coupling as defined in the Thomson 
limit. For the process at hand, this means a the electromagnetic coupling, the masses of the W, M\y, the 
Z, Mz and the top mj. 
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(a) ZZ -> tt 




(b) W~W + -> tt 



FIG. 6: Some one loop Feynman diagrams for the processes a) ZZ — > tt and b) W W + — > tt. 

The one- loop amplitudes consist of the virtual corrections •A^yoo^ an d the counterterm contri- 
butions Mct- Fig. \6\ shows a selection of the virtual corrections. There are about 500 diagrams for the 
process W~W + -> tt and 400 for the process ZZ -> tt. Though A4^J" QCD +A4aT should be ultraviolet 
finite, photon and gluon virtual exchange leads to infrared divergences. These are regulated by a small 
photon or gluon mass. For the latter the procedure is admissible as we are essentially dealing with an 
Abelian process where the triple gluon vertex for example does not show up. The photon and gluon mass 
regulator contribution contained in the virtual correction should cancel exactly against the one present 
in the photon and gluon final state radiation. A selection of diagrams in the latter category is listed in 

Fig.m 





(a) ZZ -> ttg 



(b) W~W+ -> tt-y 



FIG. 7: A selection of the final state radiation diagrams for the process a) ZZ — ► ttg and b) 
W~W + -> tt 7 . 



The photonic (gluonic) contribution is in fact split into a soft part, where the photon (gluon) energy 
is less than some small cut-off k c , Ad^g {Ej,g < k c ) and a hard part with A4y a g rd (E~ fig > k c ). The 
former requires a photon/gluon mass regulator. We use the usual universal f actor ised form with a simple 
rescaling for the case of the gluon correction. To summarise, the full one-loop electroweak correction is 
the result of all these contributions: 



\Mi\ 



2,s+v 



\M_ 



\^ EW = 2KeM B orn(M?Zp + Mar)* + \Mf l \ 2 +\M* ard \ 2 



(3.6) 



The QCD correction is defined along the same line. The percentage change, for the total electroweak 
correction including hard photon radiation will be defined as 5 EW . For the QCD correction we will define 
this as 5 QCD . 



3.2 Checks on the calculation 

i) We first check the ultraviolet (UV) finiteness of the results. This test applies to the whole set of 
the virtual one-loop diagrams. The ultraviolet finiteness test is performed by varying the dimensional 
regularisation parameter Cjjv = V e — 7£ + l°g47r, n = 4 — 2e is the dimensionality of space-time. We 
vary Cjjy by seven orders of magnitude with no change in the result. We content ourselves here and all 
other tests with double precision. The Cuv parameter could then be set to in further computations. 
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ii) The test on the infrared (IR) finiteness is performed by including both the loop and the soft 
bremsstrahlung contributions and checking that there is no dependence on the fictitious photon mass 
A 7 or gluon mass X g . 

iii) A crucial test concerns the gauge parameter independence of the results. Gauge parameter indepen- 
dence of the result is performed through a set of five gauge fixing parameters. For the latter a generalised 
non-linear gauge fixing condition [22, 25] has been chosen, 



1 

H- -J~vr i p.j ■ ■ • •»( 2' 



Cgf = - — \(d^ - ieaA^ - igc w pZ^)W^ + + i w ^{v+~5H + i~K X z)x ] 



--L(d.Z + £ z -?-(v + eH)x3) 2 -^-{d.A)\ c w = M w /M z . (3.7) 

The x represents the Goldstone. We take the 't Hooft-Feynman gauge with £jy = £z = £4 = 1 so that no 
"longitudinal" term in the gauge propagators contributes. Not only does this make the expressions much 
simpler and avoids unnecessary large cancellations, but it also avoids the need for high tensor structures in 
the loop integrals. The use of the five parameters, a, (3, S, k, i is not redundant as often these parameters 
check complementary sets of diagrams. It is important to note that in order to successfully achieve this 
test one should not include any width in the propagators. In fact our tree- level results do not include 
any width. Because of the parameters and the energies we consider, no width is required to regularise 
the cross section. We verified that the inclusion of a width, at tree-level, for the Higgs is minuscule. At 
one-loop order, effective running widths are generated. 

iv) Phase space integration for the 2^2 process is performed using the Gauss-Patterson method. While 
for integration connected to 2 —* 2 + 7 processes, we use VEGAS adaptive Monte Carlo integration 
package [26]. For the radiative process, we choose k c small enough and verify the calculation for the total 
unpolarised cross section with the help of CompHEP |18j . We then check the stability and independence 
of the result for the total and differential corrections, Eq. 13.61 with respect to k c . 

3.3 Extraction of some photonic corrections, genuine electroweak corrections 

The QED corrections contained in the electroweak corrections may be important corrections however 
they do not reveal much about the inner structure. We do not expect their contribution to be large as in 
the case of the e + e~ initiated processes which are dominated by collinear singularities. This is because 
of the large mass of the W and top. Nonetheless one can attempt to extract their effect, assuming such a 
split between weak and QED is possible. For W~W + — * ti this is not trivial because of the non-abelian 
nature of the photon radiation from the W. This said the extraction of the leading log, k c dependent, 
contribution in the virtual+soft can be considered as a genuine QED correction. For ZZ — * ti this QED 
correction can be extracted most easily by a simple adaptation of the one loop initial state radiation, 
ISR, in e + e~. We thus take, with Qt = 2/3 the charge of the top, the relative QED correction to be: 



At threshold, (3 — > 0, the QED corrections can be extremely large. This is due to the Coulomb 
correction. This can be extracted from: 
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"Coul 
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Sp 



2/3 
(3-1 



Sp 



2(3 
(3 + 1 



n 2 

+ Y 



Sp(z) 



dt 



LogQ. - t) 



(3.9) 



This Coulombic correction is also present in W W + — ► tt. The corresponding QCD correction is 
obtained by Q 2 a — ► Cfo: s , with Cf = 4/3 or in effect a — > 3a s . 
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For ZZ — ► ii the "purely" electroweak correction (PEW) is defined from the contribution of the 
soft and virtual corrections, i.e. from |A^i| 2 ' s+,; in Eg. I3.6I and therefore does not include hard photon 
radiation. To this we subtract the QED correction given by Eq. I3.8I so that: 

PEW ^a+v xqed,S+V Rom rPEW Porn / Q in\ 

&ZZ = °zz - ° °zz = ° a zz ■ W-10J 

This expression is gauge invariant, UV and IR finite and independent of the photon energy cutoff. 

For W~W + — > tt, S qed,s+V involves radiation from the initial, the final as well as their interference. 
Since we are only seeking the k c dependent term, we decided to extract it numerically. To achieve this 
we take 2 different values of the cut k c , say k c \ and k C 2- Then, for the integrated cross section the purely 
electroweak correction is defined as: 



« = a^ w (k c )-A(VS)H2k c /VS\ A(V- S ) = a ^ C i ^Y kcl > . (3.11) 

ln(fc c 2/«!clj 



3.4 Genuine electroweak corrections in the scheme 

As stated earlier we work in a scheme where the electromagnetic coupling, a, is defined in the Thomson 

limit while our processes occur at the TeV scale. The running a can induce large corrections due 

essentially to the light fermions. Moreover some ml universal correction from Ap can also be induced. It 

is therefore more advantageous to work in the G^ scheme. This helps extract the rather large universal 

corrections contained in the two point function of the vector bosons self-energies through Ar which is 

of order 3% for the Higgs masses we are considering. In this scheme the electroweak corrections are 

obtained by subtracting 2 Ar from the percentage change for the PEW correction S PEW , hence defining 
§P ew = § pew _ 2A r . 



4 One loop electroweak and QCD corrections: results 

4.1 Input parameters 

In getting our results we used the following parameters: 

a" 1 = 137.0359895 M z = 91.1875 GeV m e = 0.51099907 MeV m T = 1.777GeV m c = 1.5GeV m t = 173.7 GeV 

M w = 80A5GeV = 105.658389 MeV m u = m d = 53.8 MeV m s = 150MeV m b = 4.7 GeV 

For a s we use a running constant a s (fi) with fi = ^Jsyy. The running is evaluated at the tree-loop 
level within the MS with 5 flavours and the normalisation a s (Mz) = 0.1172. We have, for instance, 
as(^/^ww = 500 GeF) = 0.09432 and a s (^s^ = 1 TeV) = 0.08776. We also set the CKM matrix to 
unity. 

4.2 Total cross section 

The effect of the radiative corrections for Mniggs = 120 GeV is shown in figure [8] for the total cross 
section for centre of mass energies of the VV system ranging from threshold to ITeV. Results for 
Mniggs = 200, 300 GeV are listed in Table [TJ First of all, let us point to a common feature we found for 
all Higgs masses we studied. The QED corrections we extracted are relatively quite small so that there is 
little difference between the full EW corrections including hard radiation and what we call the "purely" 
EW corrections. This difference is largest, around 2% at threshold energies (around 400GeV) in the case 
of W~W + — > ti . Otherwise the difference is less than 1% for both ZZ — > ti and W~W + — > ti. Therefore 
in the remainder of our discussion we will refer to the EW corrections having in mind the full EW. For 
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FIG. 8: Total cross section for the Born, the NLO with full electroweak (EW), QCD, purely electroweak 
(PEW) with the extraction of the QED corrections and in the scheme(pew), see text, a) ZlZl — > ti 
and b) W£W^ — ► ti; Muiggs = 120 GeV. c) and d) give the percentage change. 
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both ZZ and WW, we find large positive corrections at the lower energies. As we discussed earlier, this 
is due, to a large extent, to the Coulomb correction. Apart from the low energy region, the corrections we 
find for the case of the WW and ZZ are within 3% of each other, be it for the EW or QCD corrections. 
At the threshold region though, the QCD corrections are large, especially for ZZ and increase somehow 
as the Higgs mass is increased from 120GeV to 300GeV. Past the threshold region the QCD corrections 
decrease quite rapidly turning negative with a correction of about —10% at yjsyy = ITeV. In contrast, 
past the threshold region, the electroweak corrections decrease extremely slowly and are almost constant 
with a value of around 10%. Especially at high energies, the electroweak corrections tend therefore to 
cancel out most of the QCD corrections. The electroweak corrections should therefore be taken into 
account on par with the QCD corrections. Of course, expressing the tree-level cross section in terms of 
does account for a large part of the EW corrections. In this scheme the corrections are in the range 
2 — 5%, slightly larger for ZZ with Mmggs = 120 GeV. However even in this scheme the corrections can 
be competitive with the QCD corrections as the latter become small in the energy range i/s = 500 — 700 
GeV. 



TIT" — TT/"! - +4- 

w L w L — > tt 


Mmggs = 200 GeV 


Mmggs = 300 GeV 


V s ww (GeV) 


00 (Pb) 




fiPEW 


gpew 


fiQCD 


CT (pb) 


fiEW 


fiPEW 


<^pew 


§QCD 


400 


144.33 


9.35 


11.17 


4.97 


6.44 


163.77 


10.25 


11.93 


5.42 


10.54 


500 


196.13 


8.23 


9.43 


3.23 


-4.99 


201.23 


9.41 


10.52 


4.00 


-3.97 


600 


188.49 


7.96 


8.91 


2.72 


-8.33 


189.66 


8.93 


9.82 


3.30 


-7.89 


700 


168.84 


7.79 


8.61 


2.41 


-9.87 


169.04 


8.67 


9.43 


2.92 


-9.59 


800 


148.57 


7.63 


8.36 


2.16 


-10.61 


148.52 


8.43 


9.12 


2.60 


-10.42 


900 


130.46 


7.48 


8.15 


1.95 


-10.96 


130.36 


8.22 


8.84 


2.33 


-10.81 


1000 


114.95 


7.33 


7.95 


1.75 


-11.07 


114.84 


8.01 


8.59 


2.08 


-10.93 


z L z L -> it 


Mmggs = 200 GeV 


Mmggs = 300 GeV 




00 (Pb) 


§EW 


xPEW 


gpew 


fiQCD 


0- O (pb) 


fiEW 


fiPEW 


zpew 


6 QCD 


400 


21.87 


9.96 


10.14 


3.94 


22.46 


98.10 


10.18 


10.36 


3.85 


30.87 


500 


44.24 


10.13 


10.46 


4.26 


2.93 


96.56 


10.19 


10.51 


3.99 


7.9 


600 


48.49 


10.15 


10.57 


4.37 


-4.06 


78.83 


10.21 


10.62 


4.10 


-1.42 


700 


46.89 


10.04 


10.53 


4.33 


-7.83 


65.39 


10.72 


11.20 


4.68 


-6.55 


800 


43.53 


9.84 


10.38 


4.18 


-10.17 


55.45 


10.78 


11.32 


4.80 


-9.69 


900 


39.81 


9.57 


10.16 


3.96 


-11.77 


47.85 


10.61 


11.20 


4.68 


-11.74 


1000 


36.23 


9.27 


9.91 


3.71 


-12.91 


41.87 


10.33 


10.96 


4.44 


-13.14 



TAB. 1: The Born total cross section cto and the relative corrections for the electroweak (5 EW ), "purely" 
electroweak (5 PEW ), the G^ scheme (5 pew ) and QCD (5 QCD ). 

4.3 One loop angular distributions 

We show the corrected angular distributions for an energy not far from threshold, y/syy = 400 GeV, and 
the other at the high energy of ITeV. The corresponding factors for the unpolarised case, defined 

to be the ratio of the one loop corrected differential cross section to the corresponding tree level one, 
are shown in Fig. [9j In the lower energy where the cross section is largest, see Fig. 01 the electroweak 
correction dominates over the QCD one. In the central region the two corrections are both positive and 
of very similar size. At high energy, the electroweak correction shows a very large effect in the forward 
region, while it is quite modest in the backward region and almost zero in the central region. The large 
correction in the forward region has to be put in perspective though. Indeed, in this region the differential 
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cross section is 4 orders of magnitude smaller than the in the backward region. This is just an effect due 
to hard photon radiation that migrates some of the events from the most populated region. In the central 
region the electroweak correction is almost nil. The QCD correction on the other hand is negative all 
throughout but even more so in the backward direction. This is at the origin of the relatively large and 
negative correction at the level of the integrated cross section. 




FIG. 9: The correction relative to the SM. for the unpolarised cross section WyWjj — » tt with 

Mmggs = 120 GeV, (a) Js^ = 400 GeV, (b) ^/sww = 1000 GeV . The SM Born distributions are 
dispalyed in Fig. 

5 Results at e + e level 

We now convolute the results found at the "partonic" level VV — ► tt to e + e~ — > tt + X. As discussed 
in the previous sections we impose the kinematical cut m t j > 400 GeV together with pj^ > lOGeV. We 
have considered a centre-of mass energy for the linear collider ranging from 500GeV to 3TeV. As we 
discussed previously, for ^s ee < ITeV, the EVBA approximation for the process should be taken with a 
grain of salt. We do this however in order to also comply with the analysis conducted in [6] and be able 
to compare the results for the QCD corrections we find with theirs. As in [6] we quantify the relative 
correction at the e + e~ level as: 

K vv = (a^°)/a^^l + 6 NL °. (5.12) 

Figure [TO] shows the corrections for Mmggs = IZOGeV and Mjjiggs = 300Gey as a function of the 
centre-of-mass energy, Js ee . Our results for the EW, PEW and pew corrections have a straightforward 
interpretation if one recalls the results obtained at the VV level. Indeed, we found that the EW cor- 
rections, as the VV energy is varied, is almost constant, with roughly the same correction for the WW 
and ZZ channel for the Higgs masses we studied. The correction amounts to about 10%, or to about 
5% in the G\i scheme. No wonder that this is roughly what we get after convolution. For the QCD part, 
the corrections depend on ^Js ee . For low energies, y/s ee = 500GeV, the QCD corrections for ZZ are of 
order 10% for Mniggs = 120GeV and almost twice as much for Mniggs = 300GeV. They do however drop 
quickly being sensibly equal to those found in WW, turning negative of order —5% at y / s^ = 3TeV. 
These results can again be guessed based on what we found at the VV level. We completely agree with 
the results of[6] for the QCD corrections, at least for Mniggs = 120GeV that we both consider. To 



11 







EW Z 






QCD 






EW+QCD 




M ui „, = 1 20 GeV 








, PEW Z 
















\ 






pew 











M,™,= 120 GeV 



QCD 

EW+QCD 



500 750 1000 1250 1500 1750 2000 2250 2500 2750 3000 

Vs„(GeV) 



500 750 1000 1250 1500 1750 2000 2250 2500 2750 3000 

Vs ee (GeV) 



(a) Kzz factor for e + e — > tte + e through ZlZl — + tt (b) Kww factor for e + e — > ttv^v^ through WlWl — > tt 



EW 
. _. QCD 

EW+QCD 



M Hio „, = 300 GeV 



\ 



M„ tul = 300 GeV 



QCD 

EW+QCD 



500 750 1000 1250 1500 1750 2000 2250 2500 2750 3000 

Vs„(GeV) 



500 750 1000 1250 1500 1750 2000 2250 2500 2750 3000 

Vs ee (GeV) 



(c) Kzz factor for e + e — > tte + e through ZlZl — * tt (d) Kww factor for e + e — > ttv e i/ e through WlWl — > tt 



FIG. 10: K factor representing the loop corrections to (a,c) e + e — ► tte + e through Z^Zl — > tt and 
(b,d)e + e~ — > tt^P process through WlWl —* tt. (a,b) are for Muiggs = 120 GeV and (c,d) for Mffiggs = 
300 GeV. 
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summarise, we see that for ^/s ee = 1.5TeV were these processes are most likely to be useful and where 
the EVBA can be trusted, in the WW channel the electroweak corrections are of order 10% while QCD 
correction is of order —5%. Even translating the electroweak results in the G\i scheme, the electroweak 
corrections are of order 4% and practically wash out the QCD corrections. 



6 Effects of anomalous top quark interactions 

As we argued in the introduction, the processes we are studying might be a powerful tool to reveal novel 
interactions of the top quarks to the vector bosons, in particular their longitudinal modes. Without 
picking up any specific dynamical model in particular, a simple way of studying such effects is to revert 
to the effective Lagrangian approach. The most general lowest order operators, beyond those in the 
SAi have been given in[27]. We will only pick up two such operators and see their effect in W~W + — ► ti. 
In the unitary gauge and picking up only the WWti couplings one can write: 

Zeff = 4// + 47/ = ^ig^tW+^W-u + ^-i{ia^)tW + ^W- u , (6.13) 

a\ is a scalar interaction which contribute to the S-w&ve and a2 is a magnetic interaction which contributes 
to the P-wave, see Appendix for the helicity amplitudes. A is the scale of New Physics associated with 
symmetry breaking which we take in accordance with[27j to be A = 47rt> ~ 3.1TeV. Reference [HI finds 
that ai^ of order 0.1 can be probed at a linear collider with - v /s ee =1.5TeV and an integrated luminosity 
of 200fb _1 in precisely the W fusion process we are studying. These bounds may be improved with the 
use of top polarisation [27], 

6.1 Effects in W~W+ -> it 

The dj we consider will be of the order 0.1 or lower as indicated by the bounds set at an 1.5TeV 
collider [T^]. They will contribute essentially through an interference effect with the SAi (tree- level) 
amplitude contribution. Apart from their effects on the total cross section, we wish to see how they 
affect the distributions. This not only can tell apart the contribution of the two operators but also 
differentiate between the new operators and the SAi prediction. Taking into account the SAi loop 
corrections one can ask whether the operators can mimick the loop corrections. For this purpose we set 
the value of the o« such that it reproduces the one-loop corrected SAA value for the unpolarised total 
cross section. For example, af" is the value of a\ that reproduces the full electroweak correction, 

a anomlew( a ewj = a^iEW). (6.14) 

Likewise, af d will represent the QCD correction and a e i w+qcd the total one-loop correction. We restrict 
ourselves to Mmggs = 120GeV. For example, for ^fsww = 400 GeV, we find a e ™ + i cd = 0.5 and a£' Jrqcd = 
0.15 while for yjsww = 1000 GeV, we get af™ = 0.02 an order of magnitude smaller, but perhaps too 
optimistic from the experimental measurement, at least for 1.5TeV collider, but not necessarily for a 3TeV 
machine. When we look at the angular distributions, see Fig. [TT| we see as expected, that in principle 
one can tell the different effects apart, especially for low energies. At ITeV WW centre of mass, one 
needs to get to the central region, where unfortunately the cross section is much smaller. Nonetheless, we 
clearly see that one should take into account the effect of loop corrections when looking for new effects. 

6.2 Effects in W~W + — > ti at e + e level: convolution 

We now convolute with the e/W distribution functions for the longitudinal mode and weigh the effect 
of the anomalous couplings with those of the radiative corrections as a function of the e + e~ centre of 
mass energy for e + e~ — ► tiv e iJ e . For illustration we take a\ = ai = 0.03. We look at the effect on, the 
"if -factor" K ww = a* 00 ?'**/** " 1 . 
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(a) Differential distributions for %/sww = 400 GeV (b) Differential distributions for %/sww ~ 1000 GeV 



FIG. 11: Angular distributions for the SM at tree and loop level and those including the anomalous 
couplings that give the same total cross section as the corrected total cross section, see text. The figures are 
for the unpolarised cross section. Mniggs = 120 GeV , for (a) ^sy/w = 400 GeV , (b) ^sy/w = 1000 GeV 

As expected, the effect due to the anomalous couplings increases with energy. The lesson one also 
gets is that though one can measure a discrepancy with a tree-level calculation due to the anomalous 
interaction, this is not necessarily the case compared to the use of a corrected SAi cross section. Moreover, 
"in real life" an anomalous contribution might be there but its effect is washed out or dramatically reduced 
due to the "contamination" from the QCD or the electroweak corrections, if these are not taken into 
account. This said, in "real life", other important issues must be taken care of such as the background 
and the use of a calculation that can go beyond the EVBA approximation. However, the study we 
have made here clearly shows that taking into account both the electroweak and QCD corrections is an 
important issue. 
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7 Conclusion 



There are plenty of arguments that point to the top playing a very special role in the mechanism of 
symmetry breaking. W~W + — > ti and ZZ — ► it at high energy provide a very nice testing ground for 
revealing features of symmetry breaking that may not be probed otherwise since they bring together the 
top and the Goldstone bosons. On the other hand within the SM and because of the large Yukawa 
coupling of the top, the electroweak virtual effects of the top are often large. We have performed a 
full electroweak correction to these processes together with the QCD corrections. The latter had been 
calculated recently and we agree with their findings. We find that the electroweak corrections can 
be large of order 10%, though smaller of order 5% if we use the scheme. However, even in this 
scheme the electroweak corrections can be of order the QCD corrections and combine such that the 
two corrections cancel each other for a certain range of energy. The combined effect should therefore 
be taken into consideration when looking for New Physics effects in these processes. To illustrate the 
role of the radiative corrections, both QCD and electroweak, we have also considered the effect of some 
anomalous VF-top interactions. Unless these are substantial enough to give deviations larger than 20% 
or so, one needs to take into account the radiative corrections. The latter can hide the effects or even 
increase it. Naturally these important processes can only occur as subprocesses at the LHC or the ILC. 
For the LHC, QCD top production and the QCD backgrounds are such as to prevent a study of these 
reactions. At the ILC, we have, like others, relied on the EVBA to study these reactions and therefore 
convoluted with the electron structure function describing the W content of the electron. Still in this 
context the interplay between the electroweak and QCD corrections is important, even though for the 
Higgs masses we considered the convolution for e + e~ centre of mass energy in the TeV range leads to the 
QCD and electroweak cancellling each other. Therefore taking into account the QCD corrections might 
not be enough to draw a meaningful conclusion in case a deviation is observed. Of course, to draw more 
meaningful conclusions we would need, beside a careful simulation with backgrounds, to attack the full 
radiative correction to e + e~ — > tiui>, which is a formidable task. This said the gist of the electroweak 
corrections should be contained at the level of the W~ W + — > ti. In this study we restricted ourselves to 
a light enough Higgs as allowed by the precision data. It would be interesting to review the situation in 
the case of a very heavy Higgs or a manifestation of the SM without a Higgs. 
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Appendix: Tree-level helicity amplitudes for W W + —> ti 

Let Ai = ±1 represent the transverse modes for the incoming W bosons, the corresponding longitudinal 
mode and A2 = ±1 the top quark polarisations. For the top quark the third component of the isospin 
component is I3 = \ and its charge Q = |. We also define 7? = 1/(1 — (3f ) = s/4M? and 7Tj the reduced 
propagator (7^ = sHi), where yfs stand for the center of mass energy of the WW system and 
111 -4 
IL > = 1' Uz = 7=M%' UH= s-Ml igg ; Ub = s(fi + Hi) + 4mg + 2 S PtPwcose (A 1} 

8 is the angle between the W~ and the t quark. 

The tree level helicity amplitudes for the SM are given by M X ™ X2 Xa Xa for the process W~{\\)W + (\2) — > 
i(A 3 )i(A 4 ) 

M s x ™_ Xu _ X2M = -\i\ 2 ^bPt(l ~ A 2 A)(1 - A 1 A 2 cos0)sin0 (A 2) 
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^(1 - A 2 &)(/% + Pt cos 0) + 2tt z (/ 3 (1 - A 2 ft) - 2Q^ sin 2 O w )Pw (A 4) 
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•Mq,a 1 ,-a 2 ,A2 = Al £*f ^ + A i A 2Cos0) [^(1 - \ 2 Pt)(PwPt + AiA - X 2 Pw + AiA 2 /?^ (A 7) 

M 2 i 

-2A 2 /3t cos 0) - 2A 2 7r z (I 3 (l - A 2 A) - 1Q— sin 2 9 W )(3 W (A 8) 

s J 



- /W -Ai,0,-A2,A 2 — - /v/l 0,Ai,-A2,A 2 



(A 9) 



kaSM 



- V 2 8111 1 

4 7t 



^j(PwPt + AiA + A 2 /3 W - AiA 2 /3^ + 2A 2 A cos 0) 

M 2 

+2A 2 7r z (/ 3 - 1Q— sin 2 
s 



(A 10) 
(A 11) 



KA SM — KA SM 

- /vl -Ai,0,A 2 ,A2 - • /vt Q,Ai,A 3 ,Aa 



(A 12) 



5 2 , 



o,o,-a 2 ,a 2 = "J 7w/sm 



^(1 - \2pt)(2\ 2 PwPt - /Ml - " 2/3t cos#) 

M 2 

-2vr z (/ 3 (l - A 2 /3 t ) - 2Q— ^ sin 2 %)(3 - fa) fa 



(A 13) 
(A 14) 



2 „,2 



• /x/l 0,0,A 2 ,A2 



y ((1 + - (1 - /3^)/3 w costf - 2f3 t cos 2 

M 2 

+7r H (l + - 2vr z (/ 3 - 2Q^ sin 2 W )(3 - ^)(3 W cos 



x 9 1 W 
-A 2 - 

4 7t 



(A 15) 



The anomalous helicity amplitudes including the operators in Eq. 16.131 we denote by M 



•A/Ai,-Ai,-A 2 ,A 2 = -^ / Ai,-Ai,A 2 ,A 2 =-A/ai,Ai,-A 2 ,A 2 = 



A/a 1 ,a 1 ,A 2 ,A2 = -r-PtiXzai + Aia 2 cos6») 

0,Ai,-A 2 ,A 2 = -AiA 2 -^-^/^^/?iy(l + AiA 2 COSI 
A/"-A 1 ,0,-A2,A 2 = -M),Ai,-A 2 ,A 2 



(A 16) 

(A 17) 

(A 18) 
(A 19) 
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MiMMM = ^w(hPt - X 2 Pw)sme (A 20) 

•A/"-Ai,0,A 2 ,A 2 = -M),Ai,A 2 ,A 2 ( A 21 ) 

AA ,o,-a 2 ,a 2 = -2a 2 ^^/%sin0 (A 22) 

AA ,o,a 2 ,a 2 = A 2 ^T^M 1 + + 2a 2 p w cos 9) (A 23) 

These expressions completely agree with the ones found in [15]. At high energy the anomalous part 
is given by 

A/"o,o-a 2 ,a 2 ^ -a 2 -^--^2-sin0, A/" , ,a 2 ,a 2 ^ A 2 — ^-(ai + a 2 cos6>) (A 24) 

This high energy behaviour agrees with the one given in[27]. The angular dependence of a 2 classifies 
this operator as P-wave, compared to a\ which is S'-wave[27j. 
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